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CRANK NICOLSON TYPE METHOD WITH MOVING MESH FOR BURGERS 
EQUATION 
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Abstract: A Crank-Nicolson Type method with moving mesh is constructed for numerical solution of one 
dimensional nonlinear Burgers Equation with Homogeneous Dirichlets boundary conditions. The scheme 
is derived by linearizing Crank-Nicolson finite difference scheme. Numerical solutions are obtained in 
different domains for different values of t and constant of diffusivity k. 


1. Introduction 

Moving mesh methods have important applications in a variety of physical and 
engineering areas such as solid and fluid dynamics, Combustion, heat transfer ,material 
science etc. Moving mesh method has become an important tool for computing singular 
or nearly singular problems such as interfaces, shock waves and boundary layers 
[HUANG . W (2001), MADZVAMUSE (2005)] and reaction diffusion systems in growing 
domains [BAINE(2001)]. As the domain changes shape, it necessitates the use of a 
computational mesh that moves and deforms to the new spatial configuration. The 
numerical investigation of these problems may require extremely fine meshes over a 
small portion of the physical domain to resolve the large solution variations. Successful 
implementation of this strategy can increase the accuracy of the numerical 
approximation and decrease the computational cost. 

The analysis of numerical methods using moving meshes has received limited 
attention. [Ferreira (1997)| analysed a moving mesh implicit Euler method for linear 
convection reaction diffusion problems in 1-d using an energy type analysis . It was 
shown that the method was only conditionally stable depending on the temporal and 
spatial smoothness of the moving mesh. An important issue related to the use of moving 
mesh techniques is to solve conservation laws called Geometric conservation Laws 
[THOMAS P.D(1979)]. [R. M. Fuzzeland(1990)] and J. G. Verwer gave their moving mesh 
techniques which is due to the moving finite difference method proposed [Dorfi & 
Drury (1987)]. Some practical aspects of formulations and solutions of moving mesh 
PDE’s are reviewed in [Huang(1991)]. One of the successful applications of the PDE’s 
with moving mesh is given in [C.Budd(1996)] to solve one dimensional PDE’s with blow 
up solution. [Dupont and co-workers (2002) Lie(2003)]| analysed certain moving mesh 
finite element methods for a model advection diffusion equation . [Zhijum Tan,Zhengru 
Zhang (2004)| solved physical PDE on a fixed mesh and the effect of the mesh 
movement is achieved through the grid restructuring. [J.A.Mackenzie (2007)] analysed 
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the stability and convergence of a finite difference discretization of a convection 
diffusion equation in 1D using moving mesh. 

Analytical solution and numerical solutions of Burgers equation with fixed 
domain are discussed by many authors in [Alireza,Hashemian (2008), Edris Dag Ali 
sahin (2007), E.Hopf(1950), H.Beteman(1915), J.D.Cole(1951), J. M. Burger(1948), 
Kanti Pandey and Lajja Verma(2011), N.Su.Australia(2004), Ronald E.mickens(2005), 
Sachin S.Wani,Sarita Thakar(2009), Skell & Berzins (1990)].Till now the solutions of 
Burgers equationwith moving mesh are not reported. In this paper we construct a finite 
difference scheme for one dimensional nonlinear Burgers equation with Moving Mesh. 
The equation is discretized by using Crank-Nicolson finite difference scheme. Since the 
resulting scheme is non-linear, the linearization of the scheme is achieved by defining 
an increment function. 

The paper is arranged as follows. In section 2 we present the model problem and 
its reformulation with respect to moving coordinate frame. In section 3, the moving 
mesh discretization and the linearized Crank Nicolson scheme is given. In section 4 we 
present numerical solution of Crank Nicolson Type method with moving mesh in 
different domains for different values of time t and constant of diffusivity k. Finally in 
section 5 conclusion is given. 


2. Model Problem 
Let T>0 and for each t € [0,7] ,Q, be a domain in R. We shall use the notation 


O= {(xt)eR?:x€0,,1<(,7)} 


We consider one dimensional non linear Burger’s Equation 


u, +uu, = ku (x,t)e Q; 
u = up, xeEQ,, t=0 


u =O, x€0Q,, t>0 


(1) 
Where k is a constant of diffusivity. 
The conservative form of Burgers equation is as follows 
l 2 
u, + JY —ku,| =0. (2) 


We assume that Qis bounded and uw, is sufficiently smooth. Let A, be a family of 
mappings, which at each t € [0,7] maps the points € of a reference or computational 
domain Q, in to the points of domain Q, at time t. Then for eacht € [0,7], 


A, :Q, >Q, x(€,t)=A,(€). 
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We assume that A, is a bijective and Q,=A, (Q,) is bounded. For a function 

g : Q; — R defined on the physical domain, the time derivative in the reference domain 
is 

og 


GS leer >R 


If u: Q, — R is regular enough, then by chain rule ú = = |y n% | 
t X 


The governing equation (2) in the computational domain therefore takes the 
form 


ù- iu, +54), —ku,, =0. 


On rewriting the spatial derivatives with respect to computational coordinate, we get 
; u U; 
xgit — (4-5) ug = (£) 
é 
The conservative form of above equation is 


(«eu)'—[(e-)u], = & (2) @) 


g 


4 


3. Finite Difference Scheme for Burgers Equation 
We assume that the domain Q, = |x,(t)}, x,(t)] is covered by a nonuniform mesh of 


N cells with 

x,(t)= x(t) < x,(¢)<------ < ty a(1)< xy ()=x,() 

The nonuniform moving mesh in physical space is assumed to be the image of a 
fixed uniform mesh covering the computational domain ©, = [o,1].via the mapping 


The magnitude of each physical cell will be denoted by 
h, (t)= x; (t)- xa (t)j=1, 2------ N and the midpoints of the cell are defined as 


1 , 
x, y(t) = z0) 20). j=1,2------ N. 
The location of the physical mesh points at time level f=+, and t =t, is well 
defined through the mapping x(é,t).To obtain the numerical approximation of (3) we 


require an approximation x"(E,t) of x(é,1) We will assume that x" (Et) is piecewise 


linear in space and time and hence 


x"(E,t) = x, y(t) F G z a) (xlt) + Xo (t))N for Sia SZ Š y2 and 
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x, /W=x, ylt,)¥¢-1) tt “iat, T 
Assuming this form for the mapping we have 


e= [tnk = 1 E i) F? Ses Sap (4) 





And 





ae 
TA Eai 


n+l 'n 


h o= (5t nl), Skoky (5) 


With these definitions, it is easy to see that x’ (£,t) satisfies 
n VA l, alils -h 
(x! ji = (xt) + l a ) ip oy) (6) 


To define the semidiscretization of (3), we use the notation uw; to denote the 





approximation of u(x",t,) , the solution u at (x(&;, tn), tn) and u” = (utu! ——-Uy Uy f. 


Forward and Backward divided differences in computational space are denoted by 
(D,u), = “mU (Du), a iee 
h; h; 
The central divided difference and average operator in space variables are given 
by 
U tu, 
(n), = (u; - 44) and m; = a 
Using these notations, the semidiscretization of (3) takes the form 
hyuy —3(z—% Lk = 
(gu), {[(« =) ul} = race D_)ulz} (7) 
Using Crank-Nicolson scheme for equation (3) we get 


ae e-DA + ((@-3)) 


h n+1 h n 
= — — uô 
Ge, Sn age 2 
-k pes O+-D-w0); | (8) 
AŽ 2 
In equation (8) we put 
upt = uf + vf. 


For smalll At the difference between the two consecutive solutions will be small 
and hence the quadratic terms in v; j can be neglected. And we get linear form of 


equation (8) as 
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vf [rat + rut, — 2krq| +v [C2 + 2rkq + 2rkp| 


+ vf [raptt — ruf — 2krp| 


: x n n+1 
= ut [rinit =i + 4krq + rž] +u [eD — (x2), — 4krq — Akrp| 
+uja4 [-rzj + ruj- + 4krp — rž] 
(9) 
At 1 1 


where ran Veen 





hj41 
We call scheme (9) as Crank-Nicolson Type method with moving mesh. 


4, Numerical Results 
Numerical solution of one-dimensional nonlinear Burgers equation (1) is 
obtained by Crank Nicolson Type Method (9) with moving mesh. The solutions are 
obtained for k=0.1 and k=0.01 , t=0.1,0.15,0.2.The increment in time and space are 
At = 0.001 & Ag =0.01 respectively. We assume that the computational mesh points 


are equispaced so that hj,, = h; Y j. The values of x? and x are evaluated from equation 
(4),(5) and (6). The solutions are obtained using Matlab. 


PHYSICAL DOMAIN FOR EXAMPLE 1 PHYSICAL DOMAIN FOR EXAMPLE 2 

















(a) (b) 

Fig. 1 Physical domain for Example 1 is (a) and for Example 2 is (b) 
Example1:-Consider one dimensional non linear Burgers equation (1) in the region 
bounded by x,(t)=1-e’ and x,(t)=e"* with Homogeneous Dirichlets boundary 
conditions and initial condition u(x,0)=u,(x)=sinax . The physical domain 


{[x,(t), x,(t)] , te[0,6]} is shown in fig.1(a).The solutions obtained from scheme (9) are 
shown in fig.(2) for k=0.1 and k=0.01. 
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x against u x against u 

















Fig.2 Solution of Burgers equation for example (1) with At = 0.001, Aé = 0.01 for a) 
k=0.1 b) k=0.01 
Example2:-Consider one dimensional non linear Burgers equation (1) in the region 


bounded by x,(t)=0 and x, (t) = e” with Homogeneous Dirichlets boundary conditions 


and initial condition is u,(x)=—2—. The physical domain {[x,(t), x,(t)] , te[0,1] } is 
oa 
shown in fig.1(b).The solutions obtained from scheme (9) are shown in fig.(3) for 


k=0.1 and k=0.01 . 


x against u x against u 


Pa : k=0.01 
Ppp 




















(a) (b) 
Fig.3 Solution of Burgers equation for example (2) with At = 0.001, Aé = 0.01 for a) 
k=0.1 b) k=0.01 


5. Conclusion 

A Crank Nicolson Type method with Moving Mesh is constructed for one 
dimensional nonlinear Burgers equation with Homogeneous Dirichlets boundary 
conditions. For this method discretization of the equation is obtained by Crank-Nicolson 
finite difference scheme. Numerical solutions are obtained for different domains for 
different values of t and constant of diffusivity k using Matlab. 
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